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We formally discuss the post-Minkowskian limit of /(J?)-gravity without adopting conformal trans- 
formations but developing all the calculations in the original Jordan frame. It is shown that such 
an approach gives rise, in general, together with the standard massless graviton, to massive scalar 
modes whose masses are directly related to the analytic parameters of the theory. In this sense, the 
' presence of massless gravitons only is a peculiar feature of General Relativity. This fact is never 

_^ , stressed enough and could have dramatic consequences in detection of gravitational waves. Finally 

the role of curvature stress-energy tensor of /(i?) -gravity is discussed showing that it generalizes the 
£SJ ' so called Landau-Lifshitz tensor of General Relativity. The further degrees of freedom, giving rise 

| , to the massive modes, are directly related to the structure of such a tensor. 
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O ' I. INTRODUCTION 

cr 

5h , Astrophysical observations of the last decade suggests the introduction of new ingredients in order to achieve a 
self-consistent picture of cosmos. In particular, the observation that Hubble flow is currently experiencing a speeding 
up has completely changed the approach to standard cosmology inducing to take into account theoretical approaches 
r— I ■ more general than the standard lore of General Relativity (GR). 

The simplest explanation of such a cosmic acceleration requires to include the cosmological constant in the 
Friedmann- Robertson- Walker cosmology (Concordance Model). This ingredient gives rise to a negative pressure 
contribution needed to balance the standard matter attractive interaction. Although the Concordance Model repre- 
sents the best fit model with respect to all samples of data coming from supernovae, large-scale structure and cosmic 
microwave radiation several conceptual problems come out to theoretically define and to give an explanation to 
the observed value of the cosmological constant . Furthermore, assuming the existence of both dark energy and dark 
matter, we should find out new fundamental ingredients capable of giving account to almost 95% of the total amount 
of cosmic matter-energy. 

Due to these difficulties, people have considered alternative approaches to GR that could be able to frame the 
observed late time acceleration and missing matter without introducing new ingredients. In this sense, higher order 
gravity and, in particular, fourth order gravity represent an interesting scheme which could, potentially, address 
5_i 1 the problems. 

Up to now, these theories have been investigated both at cosmological scale and in the weak field limit with 
significant results It has been shown that an accelerating late time behaviour can be easily recovered Q and, 

in addition, it can be coherently related to an early time inflationary expansion [8]. Furthermore, such an approach 
seems to deserve attention even at smaller scales. In fact, modifying the gravity action in favor of a non-linear 
Lagrangian in the Ricci scalar implies, in the Newtonian limit, corrections to the gravitational potential which can 
induce an astrophysical phenomenology interesting at galactic scales. In particular one can fit the rotation curves 
of spiral galaxies and the haloes of galactic clusters without the introduction of dark matter [9j. Besides, several of 
these extended models evade Solar System tests so they are not in conflict with positive experimental results of GR 

CM- 

A relevant aspect of higher order gravity theories is that, in the post-Minkowskian limit (i.e. small fields and no 
prescriptions on the propagation velocity), the pro pagation of the gravitational fields turns out to be characterized 
by waves with both tensorial and scalar modes [12|, LL3| • This issue represents a striking difference between GR and 
extended gravity since, in the standard Einstein scheme, only tensorial degrees of freedom are allowed. As matter of 
factSjthe gravitational waves can represent a fundamental tool to discriminate between GR and alternative gravities 

Mm 

In this paper, we want to develop, formally, the post-Minkowkian limit of analytic /(i?)-gravity models which, in 
our opinion, has never been pursued with accuracy stressing enough some peculiar points. As shown by the same 



2 



authors for the Newtonian limit, we will show that it is different from the same limit of GR since massive modes 
naturally come out in the gravitational radiation. This occurrence has a deep meaning since points out that the 
presence of massless modes only is nothing else but the particular case of GR while massive and ghost modes are 
present in general (l3j . 

The layout of the paper is the following. In Sec. HH we discuss the post-Minkowskian limit of /(i?)-gravity. 
Considerations on gravitational wave massive modes are developed in Sec lIIII Sec lIVl is devoted to the discussion of 
the role of the stress-energy tensor in /(i?)-gravity. Concluding remarks are drawn in Sec. [Vj 

II. THE POST-MINKOWSKIAN LIMIT OF f(R) - GRAVITY 

Any theory of gravity has to be discussed in the weak field limit approximation. This "prescription" is needed to 
test if the given theory is consistent with the well-established Newtonian Theory and with the Special Relativity as 
soon as the the gravitational field is weak or is almost null. Both requirements are fulfilled by GR and then they can 
be considered two possible paradigms to confront a given theory, at least in the weak field limit, with GR itself. In 
(l7l H8| . the Newtonian limit of /(i?)-gravity is investigated always remaining in the Jordan frame jlij . From our 
point of view, this is important since, by perturbatively approximating a field, some conformal features could be lost. 
Here we want to derive, formally, the post-Minkowskian limit of / (-R)-gravity. 

The post-Minkowskian limit of any theory of gravity arises when the regime of small field is considered without any 
prescription on the propagation of the field. This case has to be clearly distinguished with respect to the Newtonian 
limit which, differently, requires both the small velocity and the weak field approximations. Often, in literature, such a 
distinction is not clearly remarked and several cases of pathological analysis can be accounted. The post-Minkowskian 
limit of GR gives rise to massless gravitational waves. An analogous study can be pursued considering, instead of the 
Hilbert- Einstein Lagrangian linear in the Ricci scalar i?, a general function f(R). The only assumption that we are 
going to do is that f(R) is an analytic function. The gravitational action is then 



A 



f(R) + X£ r 



(1) 



where X = — ; — is the coupling, C m is the standard matter Lagrangian and g is the determinant of the metric. The 

c 4 

field equations, in metric formalism, read 1 

1 X 

f'Rtxv - -^fg^u - f'.^v + 9nvOgf = ( 2 ) 

dDf' + f'R-2f = yT, (3) 



with T uv = (V g£"m) energy momentum tensor of matter (T is the trace), /' = 'l^—l and □„ = . CT ;<T . We 

y/—g og^ v dR 
adopt a (+,—,—,—) signature, while the conventions for Ricci's tensor is i? M „ = R° '^^ and R a p tJ , L , = r^j; + ... for 
the Riemann tensor, where 

T aP = ^9^(9^,0 + 90a,a ~ 9<xP,o) (4) 

are the Christofiel symbols of the g^ v metric. Actually, in order to perform a post-Minkowskian limit of field equations, 
one has to perturb Eqs. (j5J| on the Minkowski background "q^u- In such a case the invariant metric element becomes 

ds 2 = g aT dx a dx T = [j] aT + h aT )dx° dx T (5) 

with small (O(h) 2 < 1). We assume that the /(i?)-Lagrangian is analytic (i.e. Taylor expandable) in term of the 
Ricci scalar, which means that 



1 All considerations are developed here in metric formalism. 
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h + .foR+^foR 2 



(6) 



The flat-Minkowski background is recovered for R = Rq ~ 0. 

Field equations ([2]), at the first order of approximation in term of the perturbation {jjij], become: 



/o 



R 



(i) 

///,' 



R^ 



-fS 



2 ^ 



(7) 



R=0 



' JO — TFT? 



dR 2 



R=0 



and □ = :a ' a that is now the standard d'Alembert operator of flat space-time. 



Where/ ° = J? . 

From the zero-order of Eqs.@, one gets /(0) = 0, while T^ v is fixed at zero-order in Eq.© since, in this perturbation 
scheme, the first order on Minkowski space has to be connected with the zero order of the standard matter energy 
momentum tensor 2 . The explicit expressions of the Ricci tensor and scalar, at the first order in the metric perturbation, 
read 



R { a} = h° 



i?« = h ar ar - Uh 



(8) 



with h — h a a . Eqs. can be written in a more suitable form by introducing the constant £ = — jj-, that is 

Jo 

- \ Uh ^ - \Kv» - ~(/w CTT - oh) v ,„ + t(afr, - v,n)(/w' CTr - nh) = ^t$> . (9) 

By choosing the transformation h^ v — h^ v — ^rj^ and the gauge condition h^" „ = 0, one obtains that field equations 
and the trace equation, respectively, read 3 



nh + 3^n 2 h = -4t(°) 

Jo 



(10) 



In order to derive the analytic solutions of Eqs. (jlOp . we can adopt a momentum- description. This approach simplifies 
the equations and allows to fix the physical properties of the problem. In such a scheme, we have : 



fc 2 M fc ) + ZfaK - k 2 Vtlu )k 2 h{k) = %T$(k) 
k 2 h(k)(l - 3£k 2 ) = jrT {0 \k) 



where 



(11) 



-ikx 



hfj,i/(k) — J (2-k)' 1 htw{ x ) 
T$(k) = f^T$(x) e-^ 



(12) 



2 This formalism descends from the theoretical setting of Newtonian mechanics which requires the appropriate scheme of approximation 
when obtained from a more general relativistic theory. This scheme coincides with a gravity theory analyzed at the first order of 
perturbation in the curved spacetime metric. 

3 The gauge transformation is h' = h^ v — — Cu,n when we perform a coordinate transformation as x' 1 * = + f M with 0(f 2 )<C 1. 
To obtain the gauge and the validity of the field equations for both perturbation hfiv and hpr, the £ M have to satisfy the harmonic 
condition = 0. 
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are the Fourier transforms of the perturbation ha V {x) and of the matter tensor T$J . We have defined, as usual, 
k x = wt — k • x and fc 2 = to 1 — k 2 ; h{k) and T^\k) are the traces of /^(fc) and T^)(k). In the momentum space, 
one can easily recognize the solutions of Eqs. (fTTj) ; h^(k) turns out to be 



r r , \ _ X T M ° (fc) £X k 2 % u kjjt^ {0) 

h ^ k >- f Q fc 2 + /^fc 2 (l-3£fc 2 ) {h { ) 

which fulfils the condition h^ v = (that is h^ u {k) fc M = 0). The perturbation variable h^(k) can be obtained by 
inverting the relation with the tilded variables. In particular, by inserting the new stress-energy tensor Sffl (fc) = 
T$(k) - i^T(°)(fc) with the trace S<°)(fc) = rfS${k), one obtains: 



X S [ °J(k) jX k 2 r]^ + 2fc A1 fc„ (0) 

fc 2 + 2.^ fc 2 (i-aefc 2 ) ( j ' ( J 

which represents a wave-like solution, in the momentum space, with a massless and a massive contributions. The 
massive term is due to the pole in the denominator of the second term: the mass is directly related with the physical 
properties of the pole itself and, thanks to the parameter £, depends on the analytic form of the model (i.e. /q and 
/q). The wavelike solution in the configuration space is obtained by the inverse Fourier transform of h^ u {k). 



III. MASSIVE MODES IN GRAVITATIONAL WAVES 



The presence of the massive term is a feature emerging from the intrinsic non-linearity of /(-R)-gravity. Specifically, 
it is related to the fact that fg ^ 0, which is zero in GR where f(R) = R. This means that massless states are nothing 
else but a particular case among the gravitational theories. A similar situation emerges also in the Newtonian limit: 
the Newton potential is recovered only as the weak field limit of GR. In general, Yukawa-like corrections, and then 
characteristic interaction lengths, are present [17| 

Some considerations are in order at this point. It is worth noticing that field Eqs. © can be written putting in 
evidence the Einstein tensor in the l.h.s. [3j. In such a case, higher than second order differential contributions can 
be considered as a sources in the r.h.s. as well as the energy-momentum tensor of standard matter: 



G,u = R,u \R9»u = + TM , (15) 

where 



rp{m) _ T, 

J- uv — 



(16) 

T (curv) _ I f(R)-f'(R)R f'(R);^-g^a g f'(R) 
J-fJ-v — 2^ f'(R) ~T f'(R) 

Actually, if we consider the perturbed metric ([5]) and develop the Einstein tensor up to the first order in perturbations, 
wc have 



G»u ~ GW - h^ u)a - -DV, - ~h,„ v - i(/W CTT - Dh)ri^ (17) 



2 ^ 2 

while the curvature stress-energy tensor gives the contributions 



TlT v) ~ e(3L " ^ v U){h ar ' aT - UK) . (18) 

This expression allows to recognize that, in the space of momenta, such a quantity will be responsible of the pole-like 
term which implies the introduction of a massive degree of freedom into the particle spectrum of gravity. In fact, 
inserting these two expressions into the the field Eqs. (1151) and considering Eqs.©, we obtain the solution: 
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nh M „(x) 



X 

'7k 



S${x) + Y,Ux) 



(19) 



where (x) is related to the curvature stress-energy tensor and is defined as 



e _ f r d 4 k k 2 r]^ + 2k fl k„ 
^ [X) ~ 2 J (2ir) 2 1 



3£fc 2 



S(°\k) 



i L-.r 



(20) 



The general solution for the metric perturbation h^ v {x) 1 when the field equations are (|15|) . can be rewritten as 



d A k S)jJ{k) , £X f d 4 k k 2 r\^ iv + 2k^k v c((n _ ; ^ , ( 



X 



fU (2^) 2 fc 2 



2/^ (2tt) 2 fc2 (1 „ 3?fc2) 



where the second pole- like term is present. In vacuum (i.e. = 0), Eqs. flD} become 

fc 2 [V(fc) + - fc 2 v)M fc )] = o 

k 2 h{k){\ - 3£fc 2 ) = 
showing that allowed solutions are of two types, i.e. : 

u = ±|k| 



= / ^^(*) e 1 ^ with fc(fc) = 



and 



(21) 



(22) 



(23) 



±Jk 2 



6 



/i(jfc) e ife:E with h(k) ^ 



(24) 



It is evident, that the first solution represents a massless graviton according to the standard prescriptions of GR while 
the second one gives a massive degree of freedom with m 2 = (3£) _1 = —^fn- Thanks to this propertiy, we can rewrite 
Eqs. fTtj]) introducing a scalar field <p = Dh so that the general system can be rearranged in the following way 



n% - - *r (0) 



3 in - 



(25) 



(□ + m 2 U = -£m 2 T(°) 
Jo 



which suggests that the higher order contributions act, in the post-Minkowskian limit, as a massive scalar field whose 
mass depends on the derivatives f'(R) and f"(R), calculated on the unperturbed background metric. 

The massive mode is directly related to the coefficients of the Taylor expansion and it is interesting to note that 
they dermine also the value of the Yukawa correction in the Newtonian approximation [l7l Il8j . On the other hand, it 
is straightforward to see that massive modes are directly related to the non-trivial structure of the trace equation as 
it is easy to see from Eq. (J3J). In GR, the Ricci scalar is univocally fixed being R = in vacuum and R oc p in presence 
of matter, where p is the matter-energy density. 
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IV. THE STRESS-ENERGY TENSOR IN /(ii)-GRAVITY AND THE GRAVITATIONAL RADIATION 

As we have seen, higher order theories of gravity introduce further degrees of fredom which can be taken into 
account by defining an additional "curvature source term" in the r.h.s. of field equations. This quantity behaves as 
an effective stress-energy tensor that can characterize the energy loss due to the gravitational radiation. Although 
the procedure to calculate the stress-energy tensor of the gravitational field in GR is often debated, one can extend 
the formalism to more general theories and obtain this quantity by varying the gravitational Lagrangian. In GR, this 
quantity is a pseudo-tensor and is tipically referred to as the Landau-Lifshitz energy- momentum tensor [20l | . 

The calculations of GR need to be extended when dealing with higher order gravity. In the case of /(i?)-gravity, 
wc have 



S / d 4 x^/^gf(R) = S / d^xCig^vjgnu^gftv.pa) 



dC 



dg P a dg pa 



pa.\ 



og P aM 



Sg 



pa 



(26) 



d 4 3 



-gH p °8g pa = 



The Euler-Lagrange equations are then 



dC 

dg7 



dC 



pa 



d 



dC 



dg P a.x ' " Xi 'dg P aM 



= 0, 



(27) 



which coincide with the field Eqs. (|2|) in vacuum. Actually, even in the case of more general theories, it is possible to 
define an energy-momentum tensor that turns out to be defined as follows : 



t A = 



oc 



dC 



dC 



gpa,a 



dg P o,\£ 



(28) 



This quantity, together with the energy-momentum tensor of matter T„„ , satisfies a conservation law as required by 

Y 

the Bianchi identities. In fact, in presence of matter, one has H pu = —T^, and then 



(V=ff£),A = -V~gH pa g P a, a = -^V~gT pa g P ^ a = -x(V=gT*), x , 



(29) 



and, as a consequence, 



(30) 



that is the conservation law given by the Bianchi identities. We can now write the expression of the energy-momentum 
tensor in term of the gravity action f(R) and its derivatives: 



dR 

dg P a,\ 



dR 



9g, 



'ptT,\£ 



dR 



9pcr,a 



9 9, 



~ 9 pa,£o 



OR 

f"R,i~ W (31) 



■dg, 



p<?M 



It is worth noticing that is a non-covariant quantity in GR while its generalization, in fourth order gravity, turns 
out to satisfy the covariance prescription of standard tensors (see also On the other hand, such an expression 
reduces to the Landau-Lifshitz energy-momentum tensor of GR as soon as f(R) = R, that is 



"OR 



dC 



9 \dg, 



GR 



pa, X 



(32) 



where the GR Lagrangian has been considered in its effective form, i.e. the symmetric part of the Ricci tensor, which 
effectively leads to the equations of motion, that is 



C 



GR 



V" 1 - 11(7 OU X LLV 1 - GO) 



(33) 
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It is important to stress that the definition of the energy-momentum tensor in GR and in /(i?)-gravity are different. 
This discrepancy is due to the presence, in the second case, of higher than second order differential terms that cannot 
be discarded by means of a boundary integration as it is done in GR. We have noticed that the effective Lagrangian 
of GR turns out to be the symmetric part of the Ricci scalar since the second order terms, present in the definition 
of R , can be removed by means of integration by parts. 

On the other hand, an analytic /(-R)-Lagrangian can be recast, at linear order, as / ~ f^R + ^(R), where the 
function J- satisfies the condition: lim^^o J~ — > R 2 ■ As a consequence, one can rewrite the explicit expression of t x 
as : 



.A _ fl.X , pi 

L a — J a |gr ' 



dR 
dg P a,\ 



--d e 



. dR 



og P a.M 



9 pa, i 



dR 
dg P a.M 



dR 

■ -r Kf- g p<Tl , 

og P a.M 



5 x a F. (34) 



The general expression of the Ricci scalar, obtained by splitting its linear (R*) and quadratic (R) parts once a 
perturbed metric © is considered, is 



(35) 



(notice that £gr = ~y/~gR)- I* 1 the case of GR £qi or ) the Landau-Lifshitz tensor presents a first non- vanishing term 
at order h 2 . A similar result can be obtained in the case of /(i?)-gravity. In fact, taking into account Eq. (f34| . one 
obtains that, at the lower order, t x reads: 



OR* 



~gpcr,£a 



— /o*a| GR + fo 



dgpa.XdJ 09paM 



f it u * dR _ l fx^u* 2 

/0%a 9pa,a 7^Ja°a rL 



g P <y,a 



(36) 



Considering the perturbed metric (0, we have R* ~ R^\ where R^ is defined as in ((8]). In terms of h and r], we get 



dR* 



()R< 



OR* 



(37) 



Clearly, the first significant term in Eq. (|3"6")l is of second order in the perturbation expansion. We can now write the 
expression of the energy-momentum tensor explicitly in term of the perturbation h; it is 



1 



t« ~ m\ aH + fo{ih^,^-nh)[h^ 4a -h^ a --S^ tpa -Dh)] 
Considering the tilded perturbation metric h pvi the more compact form 



(38) 



t = - 

"1/ 2 



-~h>\nh- -h a ah' x 

2 a 2 ,a 



(39) 



is achieved. As matter of facts, the energy-momentum tensor of the gravitational field, which expresses the energy 
transport during the propagation, has a natural generalization in the case of /(-R)-gravity. We have adopted here the 
Landau-Lifshitz definition but other approaches can be taken into account [2l| . The general definition of t x , obtained 
above, consists of a sum of a GR contribution plus a term coming from /(i?)-gravity : 



(40) 



However, as soon as f(R) = R, we obtains t x = ^q| gr - As a final remark, it is worth noticing that massive modes of 
gravitational field come out from since Oh can be considered an effective scalar field moving in a potential: t x , 
in this case, represents a transport tensor. 
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V. CONCLUDING REMARKS 



In this paper, we have formally studied the post-Minkowskian limit of /(i?)-gravity developing all the calculations 
in the Jordan frame. The main result is that, beside standard massless modes of GR, further massive modes emerge 
and they are directly determined by the analytic parameters of /(-R)-gravity, that is the coefficients /g and /g of the 
Taylor expansion. This fact is extremely relevant since it does not depend on the considered /(i?)-model but it is 
a general feature that can be enounciated in the following way: Massless gravitons are a peculiar characteristic of 
GR while extended or alternative theories have, in general, further massive or ghost states [13|. It is worth noticing 
that several indications in this sense are present in literature (l2l . I22I . [23j but their relevance, from an experimental 
viewpoint, has never been stressed enough. 

On the other hand, a similar result comes out also in the Newtonian limit of the same theories: Yukawa-like 
corrections to the gravitational potential emerge in general and they are absent only in the case of GR. It is interesting 
to note that also the characteristic lengths of such corrections are related to /g and /g' as shown in [13] • Also in this 
case, the Newtonian potential, coming from the weak field limit of GR, is only a particular case. 

These results pose interesting problems related to the validity of GR at all scales. It seems that it works very well at 
local scales (Solar System) where effects of further gravitational degrees of freedom cannot be detected. As soon as one 
is investigating larger scales, as those of galaxies, clusters of galaxies, etc., further corrections have to be introduced 
in order to explain both astrophysical large-scale dynamics 0, Q and cosmic evolution 0, [h| ■ Alternatively, huge 
amounts of dark matter and dark energy have to be invoked to explain the phenomenology, but, up today there are 
no final answer for these new constituents at fundamental level. Furthermore, the fact that, up to now, only massless 
gravitational waves have been investigated could be a shortcoming preventing the possibility to find out other forms of 
gravitational radiation. Tests in this sense could come, for example, from the stochastic background of gravitational 
waves where massive modes could play a crucial role in the cosmic background spectrum 24| 25 1 . 



[1] T. Padmanabhan, Phys. Rept. 380, 235 (2003); P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75, 559 (2003); V. Sahni, Lect. 
Notes Phys. 653, 141 (2004). 

[2] S. Nojiri and S.D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007); S. Capozziello, M. Francaviglia, Gen. Rel. 

Grav.40, 357 (2008); T. P. Sotiriou, V. Faraoni arXiv: 0805.1726 [gr-qc] (2008). 
[3] S. Capozziello, Int. J. Mod. Phys. D 11, 483, (2002); S. Capozziello, S. Carloni, A. Troisi, Rec. Res. Develop. Astron. 

Astrophys. 1, 625 (2003), arXiv:astro- ph/0303041; S. Capozziello, V.F. Cardone, S. Carloni, A. Troisi, Int. J. Mod. Phys. 

D, 12, 1969 (2003). 

[4] S. M. Carroll, V. Duvvuri, M. Trodden and M. S. Turner, Phys. Rev. D 70, 043528 (2004); S. Nojiri, S.D. Odintsov, Phys. 
Lett. B 576, 5, (2003); S. Nojiri, S.D. Odintsov, Phys. Rev. bf D 68, 12352, (2003); S. Capozziello, V. F. Cardone and 
A. Troisi, Phys. Rev. D 71, 043503 (2005); S. Carloni, P.K.S. Dunsby, S. Capozziello, A. Troisi, Class. Quant. Grav. 22, 
4839 (2005). 

[5] G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, L. Sebastiani, S. Zerbini, Phys. Rev. D 77, 046009 (2008); G. Cognola, 

S. Zerbini, Int.J.Theor.Phys. 47, 3186 (2008). 
[6] G.J. Olmo, Phys. Rev. Lett. 95, 261102 (2005); S. Capozziello, A. Troisi, Phys. Rev. D 72, 044022, (2005); S. Capozziello, 

A. Stabile, A. Troisi, Mod. Phys. Lett. A 21, 2291 (2006); T. P. Sotiriou, Gen. Rel. Grav. 38; T. Chiba, Phys. Lett. B 

575, 1 (2005). 

[7] R. Saffari, S. Rahvar, Phys. Rev. D 77, 104028 (2008). 

[8] S. Nojiri and S. D. Odintsov, Phys. Rev.D 77 , 026007 (2008); S. Nojiri, S. D. Odintsov and P. V. Tretyakov, Prog. Theor. 
Phys. Suppl. 172 81 (2008). 

[9] S. Capozziello, V.F. Cardone, A. Troisi, JCAP 0608, 001 (2006); S. Capozziello, V.F. Cardone, A. Troisi, Mon. Not. 
Roy. Astron. Soc. 375, 1423 (2007); Y. Sobouti, Astron. Astroph. 464, 921 (2007); S. Mendoza and Y.M. Rosas- Guevara, 
Astron. Astroph 472, 367 (2007); Ch. G. Boehmer, T. Harko, F. S. N. Lobo, Astropart. Phys. 29 386 (2008); S. Capozziello, 

E. De Filippis, V. Salzano, Mon. Not. Roy. Astron. Soc. 394, 947 (2009). 
[10] L. Pogosian, A. Silvestri, Phys. Rev. D 77 , 023503 (2008). 

[11] S. Capozziello, S. Tsujikawa Phys. Rev. D 77, 10750 (2008); S. Capozziello, M. De Laurentis, S. Nojiri, S.D. Odintsov, 

Gen. Rel. Grav. 41, 2313 (2009). 
[12] K.S. Stelle, Gen. Rev. Grav. 9, 343, (1978). 

[13] A. Nishizawa, A. Taruya, K. Hayama, S . Kawamura, M. Sakagami, Phys. Rev. D 79, 082002 (2009); C. Bogdanos, S. 

Capozziello, M. De Laurentis, S.Nesseris. larXiv:0911.3094l [gr-qc] (2009). 
[14] G. Magnano, L.M. Sokolowski, Phys. Rev. D 50, 5039 (1994); V. Faraoni, E. Gunzig, Int. J.Theor.Phys. 38, 217 (1999). 
[15] M. Maggiore, Phys. Rep. 331, 283 (2000). 

[16] C. M. Will, Living Rev. in Relativity 9, 3 (2006) ( [http://www.livingreviews.org/lrr-2006-3l or [gr^qc/0510072 I 



[17] S. Capozziello, A. Stabile, A. Troisi, Phys. Rev. D 76, 104019 (2007). 
[18] S. Capozziello, A. Stabile, Class. Quant. Grav. 26, 085019 (2009). 



9 



[19] S. Capozziello, A. Stabile, A. Troisi, Class. Quant. Grav. 25, 085004 (2008). 
[20] L. Landau and E.M. Lifsits, Field Theory, Pergamon Press, (1973). 

[21] T. Multamaki, A. Putaja, E. C. Vagenas and I. Vilja, Class. Quant. Grav. 25, 075017 (2008). 
[22] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). 
[23] R. Kerner, Gen. Rel. Grav. 14, 453 (1982). 

[24] S. Bellucci, S. Capozziello, M. De Laurentis , V. Faraoni, Phys. Rev. D 79, 104004 (2009). 
[25] S. Capozziello, M. De Laurentis, M. Francaviglia, Astropart. Phys. 29, 125 (2008). 



